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Abstract. We extend the standard model by adding an extra generation of isosinglet 
up- and down-type quark pair which engage in weak interactions only via mixing with 
. the three ordinary quark families. It is shown that the generalized 4x4 quark mixing 

matrix, which is necessarily nonunitary, leads to nonvanishing flavor changing neutral 
currents. We then proceed to investigate various distributions and total branching ratio 
of the inclusive B — ► X s i + £~ (I — e, fi) rare B decays in the context of this model. 
It is shown that the shapes of the differential branching ratio and forward-backward 
asymmetry distribution are very sensitive to the value of the model parameters which 
are constrained by the experimental upper bound on BR(B — > X s /i + ^). We also 
indicate that, for certain values of the dileptonic invariant mass, CP asymmetries up 
, to 10% can be obtained. 
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B factories, CLEO III and other dedicated B experiments are expected to observe 
new rare B decay channels and to improve the precision of those which have already 
been measured. Radiative B decays, which proceed via loop effects, are quite crucial 
for observing the signals of new physics beyond the standard model (SM) in the 
near future. One example of such models is the extension of the SM to contain an 
^ ■ extra generation of isosinglet quarks [1]: 

» « =1-4, (1) 



where 



u A \ _ (U 
d A )-\D 



(2) 



In other words, both left- and right-handed components of the additional pair of 
quarks, which are denoted U and D with charges +2/3 and —1/3, respectively, are 
SU(2)l singlets. As a result, the Dirac mass terms of vectorlike quarks, i.e., 
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mu(U L U R + U R U L )+m D (D L D R + D R D L ) , (3) 

are allowed by electroweak gauge symmetry. However, for ordinary quarks, gauge 
invariant Yukawa couplings to an isodoublet scalar Higgs field 0,i.e., 

- Hf'-'j.uj.o + H.C. , i = 1, 2, 3 , (4) 

are responsible for the mass generation via spontaneous symmetry breaking. The 
link between the vectorlike and the ordinary quarks are provided by the extra gauge 
invariant Yukawa couplings which can be written as follows: 

-f d ^ L D R <j> - fffaUBfi + H.C. . (5) 

After spontaneous symmetry breaking, Eqs. (4) and (5) along with Eq. (3) lead to 
4x4 mass matrices for the up- and down-type quarks: 

dlMfd R + u a L Mfu R + H.C. , (6) 

In general, Md and M u are not diagonal, and to achieve diagonalization, unitary 
transformations from weak to mass eigenstates are necessaryi.e., 

a _ au a/3 I P m _ Ad a P j/ P /y\ 
L,R L,R U L,R i u L,R — ^ i L,R U L,R ) I ' J 

where we use prime to denote the mass eigenstates. The interesting property of 
the vectorlike quark model (VQM) is that the above transformations result in 
the intergenerational mixing among quarks not only in the charged current sector 
but also in the neutral current interactions. For example, the charged current 
interaction term which, when written in terms of weak eigenstates, involves only 
the three generations of ordinary quarks, i.e., 

Jcc" = ^l^jf L ^d\W; + H.C. , (8) 

transforms to 

Jcc" = E Au'lV^d^W; + H.C. , (9) 

a,P=l V Z 

where 

= E {AlT{A d L f , (10) 
i=i 

when expressed in terms of the mass eigenstates. It is straightforward to show that 
the 4x4 quark mixing matrix V is nonunitary 



(yWf P = 5 a(S - (A d L 4a )*A d L 4p , 

{VV^f 13 = 5 aP - (A u L 4a )*A u L 4fS . (11) 

Therefore, a close examination of the neutral current sector reveals that flavor 
changing neutral currents (FCNC) like 

J Z NC " = /-V E (nU^q'hW/ - Q q sm 2 W/) , (12) 

cos v w a /3=1 

where 

jja/3 = /^<?«*\*^<?*/3 = $a(3 _ ijp^jpW = f (V^V)^ , g = down - type 
h L 1 (^ f ) Q/? , 9 = up - type 

(13) 

are consequently developed in the VQM. It is exactly for this reason that rare B 
decays can serve as excellent venues to see the effects of the vectorlike quarks [2] . 
The effective Lagrangian for B — > X s i + i~ is the following 

L cS = Q= (AsLpbeife + BsLpMRPe + 2m b CsT^b£Y^) , (14) 

where 

L„ = 7„(1 - 75) , = 7m(1 - 75) , (15) 



T» = ia^(l + l5 )q»/q 2 . (16) 

The coefficients A and B receive contributions from long-distance charm-quark loop 
(cc continuum) and the intermediate ip and ip' resonances as well as short- distance 
tree and one- loop diagrams (Fig.l). However, C, the coefficient of the magnetic 
moment operator, is purely short- disctance. The details of the contributing terms 
to the effective Lagrangian can be found in Reference [3] . The total branching ratio 
of the dileptonic B decay is dominated by the resonance contributions. However, 
by using cuts in the differential branching ratio around ip and tp' invariant mass, 
one can get information on the contributing short-distance physics [4]. 

The VQM parameters which appear in our calculations are: the U-quark mass 
m;/, the nonunitarity parameter U sh = \U sb \e ld , where 6 is a weak phase, and 

V£ s V4& = U sb — (Vc KM Vckm) • Vqkmj which is the 3x3 submatrix of the matrix 
V, consists of the elements representing mixing among three ordinary generations 

of quarks. (Vc KM Vckm) presents the deviation from unitarity of the 3-generation 
CKM mixing matrix in the VQM context. The experimental upper bound BR(B — > 
X s fi + fi~) < 5.8 x 10~ 5 [5] is used, to put a rough constraint on the magnitude 



of U sb ,i.e., \U sb \ ~ 10 3 . Since the tree level diagram proportional to U sb (Fig. 
(la)) contributes to the rare dileptonic decay channels, the one-loop terms with 

coefficient V4* V± b are significant only if Vl s V ib — (Vc KM V C km) 3> U sh . Thus, we 

parametrize our results in terms of (V^ km Vckm) /\Vcb\ — ^ instead of V^V^, 

where e = | ( Vq KM Vqkm) I/IK&I an d is another weak phase of the model. 

We use \V CS \ « 0.97, |K b | « 0.04, \V ts \/\V cb \ w 1.1, and V^Kj « [6], which 
are extracted from various experimental measurements and are not affected by the 
presence of the new physics. V* s V cb is taken to be real, as is the case, to a good 
accuracy, in the "standard" parametrization of the CKM matrix. Consequently, 
V t * s V t b, which is not known experimentally, can be expressed in terms of the VQM 
parameters as: 

\v cb \ ~ ee - lVcsl ■ (17) 

Figure 2 illustrates the differential branching ratio for some values of the VQM 
parameters as campared to the SM prediction. We use e = 0.3 and all construc- 
tive/destructive contribution possibilities of the extra, beyond the SM, terms are 
considered. We observe that away from the resonances, where SD operators are 
dominant, the shift from the SM expectation, depending on the parameter values, 
can be quite significant. To constrain the model parameters by using the experi- 
mental results on BR(B — > X s fi + fi~) reported in Ref. [5], we calculate the total 
branching but excluding the resonances ip and ip' with a 5 = ±0.1GeV cut. Our 
results are depicted in Fig. 3 in the form of acceptable regions in the \U sb \ versus 
plane for various choices of the relative sign of the extra contributions. The 
most stringent constraint is obtained if the relative phases 9 and <fi both vanish 
(Fig. 3(a)). 

Figure 4 illustrates our results for the forward-backward asymmetry distribution 
of the decay B — ► X s fi + fi~ in the VQM as compared to the SM. From Figs. 4(a) 
and 4(b) we observe that, when the tree level FCNC contributes constructively, 
even though the sign of the asymmetry remains the same as in the SM, its shape, 
away from the resonances, can be significantly different. On the other hand, as is 
shown in Figs. 4(c) and 4(d), for the destructive contribution of the nonunitarity 
induced tree level term and large enough values of \U sb \ or mjj, the sign of Aps 
can be opposite to what is predicted by the SM. One important observable in this 
decay channel is the point of zero asymmetry in the forward-backward asymmetry 
distribution, which occurs somewhere below the resonance ip. Our investigation 
reveals that, as far as the VQM is concerned, the position of this point is quite 
stable and is not shifted very much from its SM value for various choices of the 
model parameters. 

Since in the VQM, the effective Lagrangian contains terms with different CP- 
odd weak phases, as well as, LD continuum and resonance contributions which are 
sources of perturbative and nonperturbative CP-even strong phases, the direct CP 
asymmetry, unlike the SM, is expected to be non-zero. The size of this asymmetry, 



which is zero for invariant dilepton masses below 2m c threshold due to the vanishing 
strong phase, depends on the interplay of the various contributing terms. Our 
results for certain choices of the model parameters are depicted in Fig. 5. From 
Figs. 5(a) and 5(b), we observe that when the weak phase is large, smaller values 
of the nonunitarity parameter lead to significantly larger CP asymmetries sensitive 
to the U quark mass, my. In fact, for a negative relative sign of the tree level 
contribution (9 = it) and mjj = 400GeV, as is shown in Fig. 5(b), asymmetries of 
the order of 10% can be achieved. 

In conclusion, we have investigated various observables of the dileptonic rare B 
decay B — > X s £ + £~ in the presence of an extra generation of vectorlike quarks. The 
measurement of this decay channel in the near future should provide more stringent 
constraints on the model parameters. 
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Figure 1: The Feynmann diagrams which contribute to B — >■ X s £ + £~ in the VQM. (a) 
The tree level contribution (b) The box diagram (c) The penguin diagrams mediated 
by charged particles (d) The penguin diagrams mediated by neutral particles. Z and 
7 coupling to the internal quarks in (c) and (d) are proportional to U 13 " and 5 13 ", 
respectively. 




Figure 2: The differential branching ratio versus dilepton invariant mass for various values 
of the VQM parameters is compared to the SM prediction. Results with \U sb \ — 10~ 4 
and \U sb \ — 10~ 3 , for two different values of the U-quark mass, (i) m[/=200GeV, and (ii) 
mf/=400GeV, are shown by solid lines and dashed lines, respectively. 
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Figure 3: The acceptable region (shaded area) of the \U sb \ versus mjj plane for various 
values of the VQM parameters obtained by using the experimental upper bound on BR(B — > 





Figure 4: The forward-backward asymmetry distribution of B — )■ A' s /i + /i~ for various values 
of the VQM parameters is compared to the SM prediction. Results with \U sh \ — 10~ 4 
and \U sh \ — 10~ 3 , for two different values of the U-quark mass, (i) m[/=200GeV, and (ii) 
mt/=400GeV, are shown by solid lines and dashed lines, respectively. 




Figure 5: The CP asymmetry distribution of B — I for various values of the VQM 

parameters. Results with \U sh \ — 10~ 4 and \U sh \ — 10~ 3 , for two different values of the 
U-quark mass, (i) m[/=200GeV, and (ii) ni( — lOOGoY. are shown by solid lines and dashed 
lines, respectively. 



